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Abstract 

We propose a generalization of Hamiltonian mechanics, as a Hamiltonian in- 
clusion with convex dissipation function. We obtain a dynamical version of the 
approach of Mielke to quasistatic rate-independent processes. Then we show that 
a class of models of dynamical brittle damage can be formulated in this setting. 
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1 Introduction 

We are interested in the modification of the Hamiltonian formalism by adding the 
subdifferential of a convex dissipation function. In the Lagrangian formalism this can 
be traced back to Rayleigh and Kelvin (cf. Thomson and Tait [2H] or Chetayev [llj). 
For the case of autonomous Hamiltonian systems with a Rayleigh dissipation function 
added see the paper of Bloch, Krishnaprasad, Marsden and Ratiu [7J. Perturbation 
analysis of Hamiltonian systems is an old and much explored subject, which is beyond 
our scope in this paper. 

Closer to our interests is Mielke |20j theory of quasistatic rate-independent pro- 
cesses. In fact one of our purposes is to reformulate Mielke theory in a dynamical 
context. From this point of view a dissipation perturbed Hamiltonian approach seems 
the most economical. 

From the viewpoint of multivalued analysis, many generalizations of Hamiltonian 
and Lagrangian mechanics have been considered, like for example Rockafellar [25J, 
Aubin, Cellina and Nohel [6] or Clarke [12]. The problem of solving a subdifferential 
inclusion of the type (|14p for a 1-homogeneous dissipation function seems to be new. 
As a general problem this subdifferential inclusion seem to fall in the class of problems 
studied in the viabilty theory, |6j or the more recent [5], but the mathematical results 
in these papers do not apply here mainly because the dissipation is 1-homogeneous. 
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More specifically, concerning the particular form - (|42p coupled with purely Hamil- 
tonian equations (|37|) . ([39p - of this subdifferential inclusion, which is relevant for dam- 
age models in continuum mechanics, it seems that there are no mathematical results 
which could be applied to this problem as a perturbed Hamiltonian problem. We thank 
to one of the anonymous referees for pointing us to the paper [26J. From our viewpoint 
the results of this paper can be seen as leading to a solution of our problem, studied 
from the Lagrangian side, that is after reformulating it as a generalized Euler-Lagrange 
equation. Nevertheless we think that the Hamiltonian structure of this problem may 
lead to interesting discretization algorithms, maybe based on symplectic integrators, 
which are known to handle correctly the energy balance even in the discretized form. 

Outline of the paper. In section 2 we propose and begin the study of a generalized 
Hamiltonian formalism, in the form of a subdifferential inclusion using a convex dissi- 
pation function. In section 3 we show that Mielke's theory of quasistatic evolutionary 
processes is the quasistatic approximation of the formalism presented here. As an ap- 
plication, in section 4 we use the formalism for a energy of the Ambrosio-Tortorelli type 
and a 1-homogeneous dissipation function and we obtain a dynamical model of brittle 
damage which may be of interest in continuum media mechanics. 

Acknowledgements. This work has been done during two visits to LMT Cachan, 
due to the kind invitation of Olivier Allix. I want to thank him for introducing me into 
the subject of delayed damage models, as well as for many constructive discussions. 

2 Generalized Hamiltonian equations with convex dissi- 
pation 

In the Lagrangian formalism we study the evolution of a system described by a variable 
q, which satisfies the Euler-Lagrange equation associated to a function L = L(t,q,q): 



where: T represents the kinetic energy, is a smooth strictly convex smooth function 
(for example quadratic, positive definite), and £ is a potential energy or stored energy. 

In the formalism of Hamiltonian mechanics we double the variables: the system is 
described by a pair (q,p) where p has the meaning of a momentum associated with q. 
Instead of the Euler-Lagrange equation, the following system of equations is used: 




(1) 



The function L is called a Lagrangian and in many situation it has the form 




(2) 



{ 



p G D g H(t,q,p) 
q = D p H(t,q,p) 



(3) 



The function H = H(t, q,p) is called a Hamiltonian. 
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Consider for simplicity that q,p G H, where Ti is a Hilbert space with scalar product 
(■,•). The equations of Hamiltonian mechanics can be written in a compact form if we 
use the notations z = (q,p) £ H xTt, J{z) = J(q,p) = (—p,q)- 

Jz - D z H(t,z) = (4) 

In particular the Hamiltonian may take the form 

H(t,q,p)=T(p) + £(t,q) (5) 

where T represents again the kinetic energy, this time expressed as a function of p. 

In this case the two formalisms are equivalent if we take T to be the Fenchel con- 
jugate of T: 

T{p) = S up{(p,q)-f(q) : q £ ?{} 



2.1 Introducing dissipation 

Consider a "dissipation function" H(q,q), convex in the second argument, and a La- 
grangian function which is a sum of kinetic and potential energies. In the particular 
case of Rayleigh dissipation the function D has the form 

K(q,q) = ^\\q\\ 2 

where || ■ || is a norm function. Then the Euler-Lagrange equations perturbed with the 
dissipation function D are, by definition: 

D q L(t,q,q) - -DgL(t,q,q) & d 4 K(q,q) . (6) 

where the d symbol denotes the subdifferential from convex analysis. 
The Hamiltonian side of ([6]) is then 

-p € D q H(t,q,p) + dqK{q,q) m 
q = D p T(p) [t) 

This motivates us to propose the following generalization of the Hamiltonian equa- 
tions ((3D in the form of a subdifferential inclusion: 

Jz - D x H(t,z) G diH{z,z) (8) 

where d-, Tt{z, z) is the subdifferential of D with respect to i: 

di1l(z,z) = {{q,p) £HxH:Vz' = (q',p r ) £HxH (9) 

K(z,z + z') > K{z,z) + (q,q') + (p,p')} . 

We shall then be interested in the following particular case: suppose that we have 
a decomposition of the state variable q = (91,92) i n t° a non-dissipative 91 variable and 
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a dissipative q 2 variable. Then the momentum variable p decomposes as p = (p\,p 2 )- 
The Hamiltonian function H is taken as follows 

H(t,qi,q 2 ,px,p 2 ) = IC(pi) + ^(Ap 2 ,p 2 ) + £(t,qi,q 2 ) (10) 

where K, is the kinetic energy energy associated to the variable pi, A is a strictly positive 
definite symmetric operator and £ is a stored energy function. The dissipation function 
takes the form 

K(qi,q 2 ,qi,q2) = p{h) (H) 

with p a convex function. 

With these choices of functions H and D the system of equations (J7|) becomes: 

-p\ = D qi £(t,q 1 ,q 2 ) 

q\ = D p K{p!) , u , 
-P2 G D q2 £(t, qi ,q 2 ) + dp{q 2 ) 

q 2 = Ap 2 

We can see the first two equations as a Hamiltonian evolution of the variables (qi,pi) 
which has (q 2 ,p 2 ) as control parameters, coupled with a pair of evolution equations 
(the last two equations in (|12p ) for the control parameters. These last two equations 
can be see as a differential inclusion: 

-p 2 - D q2 £(t,qx,q 2 ) € d P (A- x p 2 ) . (13) 

Interesting particular cases of dissipation function p are: 

(a) p = 0, no dissipation, this corresponds to classical Hamiltonian equations, 

(b) p(q) = — (q, q), (where (•, •) is a scalar product), which can be traced back to the 
Rayleigh dissipation function, 

( c ) p(q) = \\q II) where || • || is a Banach space norm, or a more general 1-homogeneous 
convex function which, as we shall explain, is related to the approach of Mielke 
and collaborators - Mielke, Theil [22], Mielke, Theil and Levitas [23] . [20] - to 
quasistatic rate-independent evolutionary processes. 

2.2 The formalism in topological vector spaces 

We shall precisely formulate relation ([8]) for a pair of locally compact topological vector 
spaces in duality, in particular this will cover the cases of Banach or Hilbert spaces. 

X and Y are topological, locally convex, real vector spaces of dual variables x € X 
and y G Y, with the duality product (•, •) : X x Y — » R. We shall suppose that X,Y 
have topologies compatible with the duality product, that is: any continuous linear 
functional on X (resp. Y) has the form x >— > (x,y), for some y G Y (resp. y \— > (x,y), 
for some x G X). 
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In this frame we don't have scalar products, neither an equivalent of the linear 
transformation J, therefore we start by introducing natural notations which make sense 
in this generality. 

We want to study generalized Hamiltonian evolutions in the space X x Y. For a 
general element of X x Y we shall use the notation z = (x,y), or similar. 

In order to properly formulate Hamiltonian equations or inclusions we need: a 
symplectic form, a Poisson bracket and a notion of subdifferential adapted in this 
setting. These will be the most natural objects one may think about and they were 
used many times before. 

We shall use notations familiar in symplectic geometry, namely: u> for the symplectic 
form, {•, •} for the Poisson bracket, Xf for the symplectic gradient of the function 
/ : X x Y — ► R (if the linear J is available then Xf = —JDf, where Df is the 
differential of /). Instead of the usual subdifferential of a convex function F we shall 
use a "symplectic subdifferential" X F. In the usual setting in Hilbert spaces we have 
J X F = d F, where dF is the well known subdifferential from convex analysis. In this 
general setting the definition of X F is obtained from the definition of dF by replacing 
scalar products with the symplectic form. 

Remark however that in this general setting the symplectic form and Poisson bracket 
have to be understood in a weaker sense than usual, let's say on a finite dimensional 
symplectic manifold. Indeed, a symplectic form is a non-degenerated 2-form which 
is closed (we renounce to the condition of being closed); a Poisson bracket is a Lie 
bracket over a algebra of functions, with supplementary properties, while here the 
"Poisson bracket" we define sends a pair of differentiable functions from Der(X,Y) to 
a function which is not differentiable a priori. 

Let us proceed with the introduction of the necessary objects. 

Definition 2.1 The space X x Y is endowed with a symplectic form: for any z' = 
(x',y'), z" = (x",y") we define the bilinear and anti- symmetric form 

u(z',z") =(x',y»)-(x",y') . 

Der(X,Y) is the linear space of functions f : X x Y — > R which are continuously 
differentiable in each argument in the following sense: there are continuous functions 
D x f : X xY ^Y and D y f :Ix7^I such that for any (x, y) € X x Y and 

(a) for all y' € Y we have 

lim- [f(x,y + ey') - f(x,y)] = {D y f(x,y),y f ) 

6 — >U £ 

(b) for all x' € X we have 

lim- [f(x + ex',y) - f(x,y)] = (x',D x f(x,y)) 

e— >0 £ 

The symplectic gradient of f G Der(X, Y) is the function Xf : X xY —> X xY defined 
by 

Xf(x,y) = (D y f(x,y),-D x f(x,y)) 
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The Poisson bracket is the bilinear, antisymmetric form 

{•, •} : Der(X, Y) x Der(X, Y) -» R XxY 
defined by: {f,g} = u(X f ,X g ). 

Definition 2.2 Let F : X x Y — > R 6e a convex Isc function. The symplectic sub- 
differential of F is the multivalued function which sends z = (x,y) G X x Y to the 
set 

XF(z) = {z' G X x Y : Vz" elxY F(z + z v ) > F{z) + u(z',z")} 



Remark that if F G Der(X, Y) and convex then we have X F = {Xp}. Indeed, 
if we use z = Xp = (D y F(x,y), —D x F(x,y)) in the definition 12.21 of the symplectic 
differential we get 

F(z + z") > F(z) + (D y F(x,y),y") + (x" , D x F(x,y)) 

which is true due to the convexity of F. Therefore Xp(x, y) € X F(x, y). The converse 
implication, that is z' G X F(x, y) implies z = Xp(x, y), is true by standard arguments 
of convex analysis. 

We propose the following generalization of Hamiltonian evolution. 

Definition 2.3 Let H : [0, T] x X x Y -> R such that for all t G [0, T] we have 
H(t, •) <E Der(X,Y), and D : (X x Y) 2 -^lU {+00} be a "dissipation function" with 
the properties: 

(a) for any z', z" EX xY we have K(z', z") > and TZ(z' , 0) = 0, 

(b) for any z £ X x Y i/ie function 1Z(z, •) zs convex, Isc. 

Then a curve z : [0, T] — > X x Y is a solution of the evolution problem with Hamiltonian 
H and dissipation D if it is derivable for all t G [0, T] (with differential denoted by z) 
and it satisfies the subdifferential inclusion: 

z{t) - X HM (z(t)) G X (K(z(t), •)) ■ (14) 



We can give an equivalent characterization for a solution, which later will lead to a 
notion of weak solution. For any / G Der(X,Y) and any derivable curve z : [0, T] — > 
X x Y we denote by / o z : [0, T] — > R the function composition of / and z, and by 

— [/ o z\ (t) the differential of this composition. 
at 

Proposition 2.4 With the notations from definition \2.3l z is a solution of the evolu- 
tion problem if and only if for any f G Der(X, Y) and for any t G [0, T] we have: 

lZ(z{t),z(t)-X f {z(t))) > K(z(t),z(t)) + ±[foz](t)-{f,H(t,-)}(z(t)) . (15) 
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Proof. For any / G Der(X, Y) and any derivable curve z:[0,T]-»Ix7we have, 
by direct computation: 

J- t lf°z}(t) ~ {f,H(t,-)}(z(t)) = -u{z{t)-X HM (z{t)),X f {z{t))) . (16) 

Let z be a solution of the evolution problem. We choose then in (|14p z" = —Xf(z(t)) 
and use CED to get (fT5j) . 

Conversely, suppose that the curve 2 satisfies (fl~5j) . For any z" £ X x Y let us 
define / G Der(X, Y) by /(z) = u;(z, z"). It is easy to see then that Xj = —z", that 

— [/ o z] (t) = u(z(t),z") and that {/, H(t, •)} (z(t)) = ui \XH(t,-)-> z ") ■ I n conclusion 
the relation ()15|) for this choice of the function / becomes the relation (|14p for z" . □ 



It is visible that the functions / G Der(X, Y) play the role of test functions in (|15p . 
Let us consider curves / : [0, T] — > Der(X, Y), which are smooth in the sense that for 

d 

any i € [0,T] there exists —f(t,z). We suppose that the Hamiltonian H : [0,T] — > 
Der(X, Y) is such a curve. For an arbitrary t G [0, T], at each r G [0, i] we put f(r, •) in 
the relation (fT5|) and then integrate with respect to r G [0, t\. We obtain the following 
relation: 



1l(z(T),z(T)-X fiTr) (z(T))) dr > / ^(z(r),i(r)) dr + 

JO 

+ - /(0,z(0))- (17) 

dr . 



Jo . 



§- t f(r,z(r)) + {f(r,.),H(r,-)}(z(r)) 



The relation (|17|) makes sense if z is differentiable almost everywhere and 

(■r 



f 1Z{z{t),z{t)) dr < +oo (18) 
Jo 



/'[ 



§- t f(r,z(r)) + {/(r,.),ff(T,-)}(*(T)) 



dr < +oo . (19) 



io 

This is leading us to the following definition of weak solution. 

Definition 2.5 Let A be a given vector space of smooth curves f : [0, T] — > Der(X, Y) 
such that the Hamiltonian H : [0, T] — > Der(X, Y) belongs to A. Then let S(D,A) be 
the space of all curves z : [0,T]-)Ix7 which are almost everywhere derivable, such 
that Diss(z, [0, T]) < +oo and smc/i i/iai i ligj) is true /or any f £ A. 

A curve z G 5(D,„4) is a weafc solution of the evolution problem if for almost any 
t G [0, T] the inclusion {1$ is true. 

Let z G S(D,A) be a weak solution. The dissipation along this solution is by 
definition the function: 

v (t)= fu^T),^.)^))) dr . (20) 
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Proposition 2.6 Let z £ S(D,A) be a weak solution and n the associated dissipation. 
Then for any t 6 [0, T] we have 



Proof. We shall use the inclusion (|14|) . which means that for any z' € X x Y and for 
almost any t € [0, T] we have 



If we take for almost any r G [0, t] z' = — i(r) and use 7Z(z, 0) = then we get 



The desired relation is obtained by integration. □ 
2.3 The 1-homogeneous case 

Suppose that X is a Banach space and Y = X*. Then X x Y is a Banach space and 
the natural norm onlxY induces a distance d(z' , z v ) = \\z' — z" \ \. 

Suppose moreover that for any z S X x Y the dissipation function D has the 
property that 1Z(z, •) is positively one-homogeneous. Then the dissipation function can 
be seen as a dissipation metric in the sense that it induces: 

(a) a "dissipation length" defined for any curve z : [0, T] — > X x Y which is almost 
everywhere differentiable by: 



The space of curves with finite dissipation length is denoted with (X x Y). 

(b) a "dissipation distance" D : (X x Y) 2 -> R U {+oo}, where D{z',z") is defined 
as the infimum of the dissipation lengths of all curves joining z' and z" . 

(c) a "dissipation variation" defined for any curve z : [0, T] — > X x Y as: 



BVji(X x Y) denotes the space of curves with bounded dissipation variation. 

The dissipation distance D is not really a distance, because it is not symmetric and 
it may take the value +oo. It satisfies nevertheless the triangle inequality. 

The dissipation length and dissipation variation are defined in principle for different 
classes of curves, but in particular cases they are the same. All in all this is a general- 
ization of well-known facts in the analysis in metric spaces, see for the relevant results 




n{z{t),z{t) + z') > K(z(t),z(t)) +u(z(t)-X H{t> . ) (z(t)),J) 



W (i(T),A^ (Ti .)(z(r))) > K(z(t),z(t)) > 
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Gromov chapter 3 [19], or Ambrosio, Gigli, Savare chapter 1[3], which has been devel- 
oped by Mielke and collaborators in the theory of rate-independent evolution systems 
(see section [3] for further details and references). Enough is to mention that if z is a 
curve which is differentiate almost everywhere and of finite dissipation length then its 
dissipation length equals the dissipation variation. 

In particular then any weak solution satisfies (|17p with the term / 1Z(z(t), z(t)) dr 

Jo 

replaced by Diss(z, [0, t]). If the class A is sufficiently rich then satisfaction of (]17p will 
imply that z is a weak solution. 

Theorem 2.7 If TL{z' , z") = TZ-(z',x") for any z',z" € X x Y then for any weak 
solution z:[0,T]->Ix7 and for any t € [0, T] we have: 

ff(0,*(0)) + J —H(t,z(t)) dr = H(t,z(t)) + Diss(z,[0,t\) (21) 



Proof. In relation (|1T|) let us take / = XH for an arbitrary A € (— oo, 1): 

f K(z(t),z(t) - \X H(iy) (z(T))) dr > Diss(z,[0,t}) + XH{t,z{t)) - XH(0,z(0)) 
Jo 



X 



d_ 
~dt 



H(t,z(t)) 



dr 



(22) 



In the hypothesis of the theorem if z is a weak solution then it satisfies the following: 
for almost any t G [0,T] and for any z" = (x",y") elxY 

n(z(t),x(t) + x") > K(z(t),x(t)) + (x(t)-D y H(t,-)(x(t),y(t)),y"}- 

-(x",y(t)+D x H)t,-)(x(t),y(t))} 

It follows that for almost any t 6 [0, T] we have x(t) = DyH(t,-)(x(t),y(t)), therefore 
for almost any r G [0, t] we have: 

7e(z(r),i(r)-AX H ( T ,.)^(r))) = 7Z.(«(r), x(r) - XD y H(r, -)(x(r), y(r))) = 

= rc(*(r),(l-A)(£(r))) = (1 - \)K(z(t), x(r)) 



We return to (|22p . we use the information that we gained and the equality between 
dissipation variation and dissipation distance and we obtain: for any A € (— oo, 1) we 
have: 



> A 



Diss(z,[0,t]) + H(t,z{t)) - H(0,z(0)) - [ 

Jo 



d_ 

dt 



H(t,z(t)) 



dr 



The arbitrary A can have any sign, therefore we deduce the desired equality (|21|) from 
the previous inequality. □ 
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This theorem shows a great advantage of Hamiltonian formulations upon Lagrangian 
formulations: a weak Hamiltonian formulation naturally conserves quantities of inter- 
est, like the energy, while in Lagrangian formulations this has to be imposed by hand 
(which then leads to different weak and energetic formulation). This can be stated in 
few words as: weak solutions of the Hamiltonian formulation are energetic solutions in 
the Lagrangian formulation. 

3 Connection with Mielke's theory of quasistatic evolu- 
tionary processes 

Consider a physical system with the state space Q. This space may have a manifold 
structure, or it may be a space of functions q : Q —> Ai, with given regularity, where 
M is a manifold. In this case the bounded Lipschitz domain 0, represents the reference 
configuration of a continuous body. We shall denote a generic point of Q by the letter 
q and q denotes a vector in the tangent space to Q at q G Q. 

For the first time in the proceedings paper Mielke, Theil [22], then in Mielke, Theil 
and Levitas [23] , the notion of a energetic solution of a quasistatic evolutionary process 
was introduced, based on a energy function 

£:[0,T]xQ^Eu{+oo} , £ = £(t,q) 

and a "dissipation metric" 

K : TQ -> [0,+oo] , K = K{q,q) 

Here TQ = {(q,q) \ q G T q Q} is the tangent space space to Q at q G Q, in a generalized 
sense. 

The dissipation metric is convex and lower semicontinuous with respect to the 
second variable. For the case of rate-independent processes the dissipation metric is 
1-homogeneous (i.e. it can really be interpreted as a metric). The force balance 
equation is: 

G dgTZ(q,q) + D q £{t,q) (23) 
To the dissipation metric 1Z is associated a non symmetric dissipation distance 

D:QxQ^ [0,+oo] 

D( qi ,q 2 ) = mf{j\(q(s),q(s))ds | q G W rl ' 1 ([0, 1], Q) q(0) = qi, q(l) = q<A 

Definition 3.1 A evolution q : [0, T] — > Q is an energetic solution associated with £ 
and D if 

(a) the functiont G [0,T] i— > d t £(t,q(t)) belongs to L 1 ((0,T)) ; and for every t G [0,T] 
we have £(t,q(t)) < +oo, 

(b) the stability condition holds: for any q G Q 

£(t,q(t)) < £(t,q) + D(q(t),q) 
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(c) the energy balance holds: 

£(t,q(t)) + Diss(q,[0,t}) = 5(0, g(0)) + f d t £(s,q(s)) ds 

Jo 

where 

f N \ 

Diss (q, [0, t]) = sup < D(q(sj-i), q(sj)) \ all partitions of [0, t] > 



We can recover the force balance equation (|23p from the generalized Hamiltonian 
formalism with dissipation proposed in section [2j Indeed, suppose that the state space 
of the physical system is Q = B, a reflexive Banach space. Consider the phase space 
X x Y = B x B*. A generic element of z £ has the form z = with g € B, 

p € B*. 

We shall take Hamiltonian and dissipation functions almost as in (|10p . (jlip . The 
Hamiltonian function has the form H(t,q,p) = fC(p) + S(t,q) where /C is a smooth 
function (kinetic energy) and £ is the energy function of Mielke. We take a dissipation 
function 7Z(q,p, q,p) = TZ(q, q) with 1Z the dissipation metric. 

With these choices of functions H and D the equation (I14h takes the form: 

-p G D q £(t,q) + dgK(q,q) 
q = D p K{p) 

The quasistatic version of (f23l) is just the force balance equation of Mielke (f23l) . We 
are also in the hypothesis of theorem I2.71 If we neglect the inertial terms in (|2ip we 
obtain the energy balance condition (c) from the definition of energetic solution I3.ll 

Let us see what is the expression of the dissipation along a solution of (plj) . as 
defined by (JUJ). We have 



m = -{D p JC(p(t)),p(t)) - {q(t),D q £(t,q(t))) 
As in the proof of proposition I2.6I we arrive to the inequality 

0> K(q(t),q(t)) + (D p )C(p(t)),p(t)) + (q(t),D q £(t,q(t))) 
therefore we get TZ(q(t),q(t)) < fj(t). We integrate this inequality and we obtain: 

> I n(q(s),q(s))ds 



o 



We finally obtain that r](t) > D(q(0),q(t)) > 0, which means that the dissipation rj 
along a solution of (|24|) is always greater or equal to the dissipation distance (in fact 
greater than the dissipation length). 
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4 Application: a dynamical model of brittle damage using 
the Ambrosio-Tortorelli functional 



Mielke and Roubfcek [21] proposed a rate-independent brittle damage model based 
on the theory of rate- independent evolutionary processes |20j . The model of Mielke 
and Roubfcek is a quasistatic particular case of the more general dynamical model of 
Stumpf and Hackl [27] . 

By using the generalized Hamiltonian formalism we are able to obtain a dynamical 
model of brittle damage, which is also a particular case of the general dynamical model 
of Stumpf and Hackl. 

The model is based on a energy of Ambrosio-Tortorelli type and a dissipation func- 
tion as in the model of Mielke and Roubfcek. 

4.1 The Ambrosio-Tortorelli functional 

Let n € N be a strictly positive natural number and f2 C M n a bounded, open set, with 
piecewise smooth boundary. The Mumford-Shah functional [21] is 



defined over all pairs (u,S) such that u G C (Q\S,R). The set S is a n — 1-dimensional 
surface in ]R n , or a countable union of such surfaces. In the case n = 2 this functional 
can be seen as the energy of a brittle body suffering a antiplane displacement u and 
presenting a crack S. 

For n = 3 the state of a brittle body is described by a pair displacement-crack, 
(u, S) is such a pair if S is a crack (a 2D surface) which appears in the body and 
u e c\n \ s,r 3 ) is a displacement of the broken body, that is u is smooth in the 
exterior of the surface S, but it may have jumps over S. The total energy of a brittle 
body is a Mumford-Shah functional of the form: 



The first term of the functional E represents the elastic energy of the body with the 
displacement u. The second term represents the energy consumed to produce the crack 
S in the body. Here his energy is taken to be proportional with the area of the crack S 
(technically this is the 2 dimensional Hausdorff measure of S), with the proportionality 
factor G, which is the Griffith constant. 

Starting with the foundational papers of Mumford, Shah [24], De Giorgi, Ambrosio 
[13] . Ambrosio [T], [2], the development of models of quasistatic brittle fracture based 
on Mumford-Shah functionals continues with Francfort, Marigo [16], [T7|, Mielke [20] . 
Dal Maso, Francfort, Toader, [13], Buliga [8], [9], [TO]. 

All these models are based on a technique of time discretization followed by a 
sequence of incremental minimization problems. These models are either seen as appli- 
cations of De Giorgi method of energy minimizing movements, or in the frame of the 
theory of Mielke of rate- independent evolutionary processes [20J . 




(25) 




(26) 
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The functional 



E c (u,d) = [<j)(d)^K | Vu | 2 + ]pc | Vd | 2 + ^ 2 j (27) 

was introduced by Ambrosio and Tortorelli [4], as a variational approximation of the 
Mumford-Shah functional (|25|) . Here d is a field which approximates the characteristic 
function of a crack, that is d : 0, — > [0, 1] and the set 

S c = {x e O : 1 > dc(a;) > 1 - 0(c)} 

approximates the crack. More precisely, if (u c , <i c ) is a minimizer of the Ambrosio- 
Tortorelli functional (|27p then as c — > the displacement ti c converges (in some norm) 
to a displacement u, the set if? c shrinks to a surface S" and (u, S 1 ) is a minimizer of the 
Mumford-Shah functional (|25p . 

The variable <i plays the role of a brittle damage variable, because it takes values 
in [0, 1] and also because it is coupled with the antiplane displacement u through the 
term 

' |2 



which represents the elastic energy of the body with elasticity coefficient <p(d)K. The 
function (ft is taken as a decreasing function from [0,1] to [0,1], such that 0(1) = 0, 
0(0) = 1. 

Focardi [15] proved that there is a Ambrosio- Tortorelli functional suitable for ap- 
proximating the 3D Mumford-Shah functional (|26|) . namely: 



E c (u,d) = jU(d)w(Vu) + ~7c | Vd | 2 +^ 2 j (28) 
under certain growth conditions on the elastic energy function w. 



4.2 Quasistatic model, using Mielke's theory 

In this subsection we obtain an interpretation of a mathematical result of Giacomini 
|18j . which shows that models of damage based on the Ambrosio- Tortorelli functional 
have the important property of being compatible with brittle damage from the energetic 
point of view. This is a desirable feature of a model of brittle damage, as there are 
many "classical" models of brittle damage which allow the creation of a brittle crack 
(seen a concentrated total damaged region) with zero consumed energy. 

We shall look at the equations coming from the force balance equation of Mielke 
(|23p and the Ambrosio- Tortorelli functional taken as the potential energy. The state 
of the system is described by a pair (u, d), where u is the displacement and d a scalar 
damage variable taking values in [0, 1]. 

We shall take a dissipation metric which is almost the same as in Mielke and 
Roubi'cek model [21], relation (2.5) (see also the discussion at the end of the section 
2.2), which gives the dissipation functional 

K(u, d, u,d) = J [f5d + xi(d) + X2(d)} 
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The functions Xl>X2 are indicator functions of convex sets: 



Xl (d) 



{ 





+00 



if d e [0, 1] 



else 



Xi(d) 





+00 



else 



if d <E [0, +00) 



Formally integrating by parts the force balance equation of Mielke (|23p . we arrive 
to the evolution equations: 



is greater or equal than 0, due to the fact that (ft is decreasing, thus eft' < 0. This term 
represents the variation of the elastic energy density due to damage. 

The paper [18] can be seen as an investigation o f the limit to the fracture model of 
the bulk damage model of Mielke and Roubi'cek, that is in the limit when the damage 
variable equals almost everywhere (therefore the value of the parameter (5 is not 
important in the sense that a > makes the same effect as /3 = 0). This result can be 
described as follows: for any parameter c let q c = (u c , d c ) denote an energetic solution 
associated with the Ambrosio-Tortorelli energy E c and dissipation distance D coming 
from the dissipation metric 7Z. Then as c converges to 0, the evolution q c converges to 
an evolution (u, S) of the energetic formulation of brittle fracture of Francfort, Marigo 
[E] or Buliga [9]. 

From the point of view of mechanics fracture is a manifestation of concentrated 
damage. Therefore a good (bulk) damage model should have the property that it is not 
possible to produce arbitrarily concentrated damage with arbitrarily small expense of 
energy. Such models are said to be compatible with brittle fracture from the viewpoint 
of energy balance. There are many models of brittle damage in use, not all of them 
compatible with brittle damage. The mathematical result of Giacomini means that the 
Ambrosio-Tortorelli functional leads to brittle damage models which are compatible 
with brittle fracture from the point of view of energy balance. 

4.3 Hamiltonian brittle damage 

We shall apply the generalized Hamiltonian approach to a functional of the Ambrosio- 
Tortorelli type. 

We take as state Q = (u, d) the pair formed by the displacement u and the scalar 
damage variable d E [0, 1]. The space of this pairs corresponds to the space X from 
the general model. 

The dual variable, in the sense of Hamiltonian mechanics, is P = (p, y) £ Y, where 
p is the momentum and y is a scalar variable dual to d (which will turn out to be 
linearly dependent on d). 




G - jcAd + <ft'{d)\K I Vu | 2 +£d + d X 2{d) 



(29) 



We add the constraints d € [0, 1], boundary and initial conditions. The term 



(ft'(d) -K I Vu 
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The space of all pairs (Q, P) is a product of two symplectic vector spaces B = 
B\ x 82- The space of non-dissipative variables B\ x B*[ is a space of pairs of (weak) 
functions (u, p) defined over fi C M 3 . Therefore u G B\ and p € B^, where B\ is a 
Banach space (for example a well chosen Sobolev space of functions over Q) and B\ is 
its dual. The duality product is 

(P,u)i = / p u 
Jn 

Similarly, the space of dissipative variables (d,y) is B2 = B2 x B%, a space of pairs of 
(weak) functions (d, y) defined over !1C1 3 , with B 2 another Banach space of functions 
over d, B\ is its dual. The duality product is 

(y,d) 2 = yd 
Jn 

Let us define the the Hamiltonian as: 



H(t,u,p,d,y) = £(u,d) + T(p,y) - (Z(i),u) (30) 
where ^ is the free energy, T the kinetic energy and l(t) the external forces, seen as: 



(l(t) 



u) = [ f(t) • u + ( f(i) • u 



Here T C <9S7 is that part of the boundary where surface forces f(t) are imposed at the 
moment t. The stored energy £ is therefore: 

£{t,u,d) = *(u,cQ - (l(t),u) . 

Displacements may be imposed on another part V of the boundary <9$7. This is done by 
imposing that at every moment t E [0, T] the displacement u(i) belongs to a subspace 
B\ (t) C B\ of kinematically admissible displacements. 
The expression of the free energy is the following: 



*(u,d) 



<f>{d)w{Vu) + ^||Vd|| 2 + ^L\d\ 2 



(31) 



which has a form analoguous with the one proposed by Stumpf and Hackl [27] formula 
(3.34). Here <j) is a smooth, decreasing function with values in the interval [0, 1]. 
The kinetic energy has the form: 



T(p,y) 



1,1 |2 1 „ l|2 

2*1 y I 2 +^IIpII 2 



(32) 



The second term in the expression of the kinetic energy is just the usual kinetic energy 
expresses as a function of momentum p, as it is usual in the Hamiltonian formalism. 
Similarly, y is a momentum variable corresponding to d and b is the scalar version of 
a microinertia tensor (we use the same name as Stumpf and Hackl [27] concerning the 
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kinetic energy described in their formula (2.4)). We suppose that the constants K, L 
and b are positive. 

The dissipation function is the same as in the previous section: 



X[o,i](d) + X[o,+oo)(d) + /3\d 



(33) 



We shall find the equations satisfied by any curve of evolution (u, p, d, y) : [0, T] — > B 
which is a solution of the generalized Hamiltonian equations ()14|) . for the Hamiltonian 
([HO]) and dissipation (f3"3"j) . By using the expressions of the free energy (f3"T|) and kinetic 
energy ([32]) . we obtain: 



(p, u)i + (Ai$(u, d), u)i = {1(f), u)i Vu G B a (t) 



(34) 



(p, u)i - (p, DpK(p, y))i = Vp G 



(35) 



There are two more equations, for the evolution of d and y. Due to the non smooth 
dissipation, these are in fact expressed as subdifferential inequalities: for almost any 
t G [0, T] d(t, x) G [0, 1] and y(t, x) G [0, +oo) for almost every x G ft, at any t G [0, T] 
the displacement u(i) is kinematically admissible, i.e. u(i) G B\(t), and moreover for 
any d G B2, such that d(x) + d > for almost every 1 G O, and for any y G -B2 we 
have: 



/3 



> 



(36) 



> 



y(d - 6y) - (Ld + ^'(d)w(Vu) + y)d - KVdVd 



The equation (|34|) gives the usual momentum balance: for any u kinematically 
admissible we have 



-p • u - 0(d) Di«(Vu) : Vu] = / f(t) • u + / f(t) • u 
n Jn Jr 

Denote by S = Dw(Vu) the stress variable given by the elastic energy w. Integration 
by parts leads us to a balance equation and boundary conditions: 

div {<P{d) S) + f(i) = p in Q (37) 

0(d) Sn = f(i) on r , 0(d) Sn = on dtt \ (T U V) , u = u (t) on V' . (38) 
Equation (|35|) gives us the momentum p as function of u: 

p = pu . (39) 

Equation (|36p is equivalent to the following two relations: 

d = by (40) 
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and for all d € B2, such that d(x) + d > for almost every x € Q 



P 



d+d 



d 



> 



(41) 



> - {Ld + (j)'(d)w(yu) + y)d + KVdVd 
Jn L 

Let S 1 : [0, +00) — > 2 M be the multivalued function defined by: 

cr \_j P , v >0 

b W-\ (-co,/?] , « = 

The function S 1 is the sub differential of a convex function. By using the definition of 
5 1 and relation ([40 p we obtain the following equivalent form of the inequality (|41|) : for 
almost every iGfiwe have: 



(y + + 0'(d)to(Vu) - iTAd) € S(y) 



(42) 



We may add the boundary condition (which is not strictly speaking a consequence of 
the formalism): on d£l we have: 



y>0, -K^deS(y) 
an 



(43) 



In the particular case of a functional of the Ambrosio-Tortorelli type (j28[) we may 
take: 

K = 7c , L = — , b = 7c 



The function which enters in the expression of the free energy is chosen as in the 
Ambrosio-Tortorelli functional. With this choice of constants we obtain from (|42|) and 
(|4U|) the differential inclusion: 

- (d + 7 2 d + jccj)'{d)w{Vu) - 7 2 c 2 A(i) G jcS(d) 

This inclusion suggests that in this model there is a maximal speed of propagation of 
damage of order 7V 1 + c 2 . 
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